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ABSTRACT 

The growth of supermassive black holes (BHs) located at the centers of their host galaxies comes mainly from 
accretion of gas, but how to fuel them remains an outstanding unsolved problem in quasar evolution. This issue can 
be elucidated by quantifying the radiative efficiency parameter (77) as a function of redshift, which also provides 
constraints on the average spin of the BHs and its possible evolution with time. We derive a formalism to link r\ 
with the luminosity density, BH mass density, and duty cycle of quasars, quantities we can estimate from existing 
quasar and galaxy survey data. We find that rj has a strong cosmological evolution: at z « 2, rj « 0.3, and by z ~ 
it has decreased by an order of magnitude, to rj sa 0.03. We interpret this trend as evolution in BH spin, and we 
appeal to episodic, random accretion as the mechanism for reducing the spin. The observation that the fraction 
of radio-loud quasars decreases with increasing redshift is inconsistent with the popular notion that BH spin is 
a critical factor for generating strong radio jets. In agreement with previous studies, we show that the derived 
history of BH accretion closely follows the cosmic history of star formation, consistent with other evidence that 
BHs and their host galaxies coevolve. 

Subject headings: black hole physics — galaxies: evolution — quasars: general 



1. INTRODUCTION 

The empirical correlation between black hole (BH) mass and 
bulge luminosity (Kormendy & Richstone 1995; Magorrian et 
al. 1998; Marconi & Hunt 2003) and stellar velocity dispersion 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000; Tremaine et 
al. 2002) have stimulated the notion that BHs evolve jointly 
with their host galaxies, mediated perhaps through feedback 
from active galactic nuclei (AGNs; e.g., Di Matteo et al. 2005). 
Many issues affect our detailed understanding of the relation- 
ship between BHs and their hosts, including (1) the evolution 
of BH mass, accretion rate, and spin, (2) the assembly and evo- 
lution of galaxies, and (3) the mutual interaction between BHs 
and galaxies. A simple comparison of the mass density of BHs 
in local galaxies with the integrated energy density of quasars 
(Soltan 1982) leads to the conclusion that radiatively efficient 
accretion, with an average efficiency of rj w 0.1, is the main 
driver of BH growth (Small & Blandford 1992; Chokshi & 
Turner 1992; Yu & Tremaine 2002; Marconi et al. 2004). If the 
angular momentum of the accreting material is aligned with the 
spin of the BH and remains constant, we expect the BH to spin 
up to its maximum rate after accreting one-third of its mass, 
and attain rj sa 0.4 (Thorne 1974), in apparent conflict with the 
efficiency deduced from Soltan's argument. This inconsistency 
can be plausibly resolved if BHs are fed by episodic, random 
accretion (King & Pringle 2006), which has a tendency to spin 
down the BH over time (King et al. 2008). If this, indeed, is the 
main mechanism for BH growth, we predict that the radiative 
efficiency of quasars should decrease from high to low redshift. 
The main goal of this paper is to test this prediction observa- 
tionally. By comparing the cosmic history of BH accretion with 
the better-established cosmic history of the star formation rate, 



we can further constrain the purported joint evolution of central 
BHs and their host galaxies (e.g., Silk & Rees 1998; Croton et 
al. 2006; Wang et al. 2007). 

In this Letter, we set up a formalism to describe the depen- 
dence of rj on the luminosity density, growth rate, and duty 
cycle of quasars. We apply it to survey data and show that rj 
evolves strongly with redshift, which we interpret as a signa- 
ture of BH spin-down as a consequence of episodic, random 
accretion. Our calculations assume a standard cosmology with 
ff = 71kms _1 MpcT 1 , Q m = 0.3, andfi A = 0.7. 

2. THE ?/— EQUATION 

Following Soltan's argument, we assume that the growth of 
BHs comes mainly from accretion during their active (AGN) 
phases. From its definition (Marconi et al. 2004), the radiative 
efficiency is 

where Ae and Ap. are the cumulative increases of the energy 
and the mass density of BHs from z+ Az to z, and c is the speed 
of light. We can define the duty cycle as the ratio of the number 
density of quasars (Nq) to that of active and inactive galaxies 
(N a ), 

N Q (z) 



S(z) = 



N Q (z)+N G (zy 



which can be equivalently expressed as (Wang et al. 
2008) 

~ qso (z) 



S(z) = 



Pi 



Pacc(z) 
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where pT°(z) is the mass density of quasar BHs and p acc (z) 
is the mass density of all gas accreted onto BHs until z. The 
mass density of quasar BHs can be obtained from pT°(z) = 
f M , M.$(z,M.)dM„ where <I>(z,M.) is the mass function of 

BHs with the break given by M* . Since quasars can increase 
their masses only via accretion (i.e. Ap. = Ap acc ), we have, us- 
ing the equivalent formulation of the duty cycle, the cumulative 
increase of the mass density 

A P . = A(<rVD- (4) 

The cumulative energy density can be obtained through 

Ae = I/(z)Pj0 Az, (5) 

where t is cosmic time, U(z) = J L TLboi^QfcLboi)dLboi is the 

bolometric luminosity density, 5 , Q(z,Lboi) is the bolometric lu- 
minosity function (LF) of quasars (Hopkins et al. 2007), and 
the lower limit Lq, given by the break in the LF, is determined 
from quasar surveys. We then have the r/-equation: 

Equation (6) is striking in that it describes 77 in terms of 
pT°/S, which combines the growth and the episodic cycles 
of BHs with the evolution of their host galaxies. This equa- 
tion has the advantage of only depending on observables: (1) 
pT°(z), which can be reasonably estimated from empirical rela- 
tions calibrated against reverberation mapping; (2) S(z), which 
in principle can be estimated from deep galaxy surveys; and (3) 
U (z), which can be directly obtained from quasar LFs. Further- 
more, the equation directly gives r\ at any redshift, whereas the 
usual Soltan's argument requires a comparison of the accreted 
mass density of BHs with the local mass density of BHs to get 
the averaged radiative efficiency over time (Soltan 1982; Yu & 
Tremaine 2002; Marconi et al. 2004). These features allow us 
to deduce the radiative efficiency as a function of redshift from 
observables. 

3. EVOLUTION OF THE RADIATIVE EFFICIENCY 

Equation (6) can be applied to existing survey data. We use 
the duty cycle expressed by Nq(z) and Nc(z), where Nq(z) = 
f L « *Q(z,L)dL and N G (z) = f M , ^ G (z,M R )dM R are the number 

densities of quasars and galaxies from their luminosity func- 
tions "Pq and respectively. It should be noted that is 
wavelength-dependent (Wolf et al. 2003), being less sensitive 
to stellar population changes at longer wavelengths (Ilbert et al. 
2005). We calculate Ng(z) from the rest-frame /?-band LFs of 
galaxies up to z ~ 2 derived from the deep HST study of Dahlen 
et al. (2005) and the VIMOS-VLT Deep Survey (VVDS) of Il- 
bert et al. (2005). Both samples have similar selection criteria 
and photometric redshift estimates from z ~ to 2. Nq(z) is 
calculated from Bongiorno et al. (2007), who supplement the 
Sloan Digital Sky Survey (SDSS) quasar sample of Richards et 
al. (2006) with VVDS data for low-luminosity AGNs. We cal- 
culate pi™ with BH masses estimated from the "virial" method, 
as given in Vestergaard et al. (2008), but use the corrected value 
for z = 0.49 and the new value for z = 0. 17 in Kelly et al. (2009). 
The samples of quasars and galaxies used in these LFs overlap 
in at least some survey fields, minimizing potential observa- 
tional biases. We relate the lower limits of the galaxy LF (M^) 
and the BH mass function (M*) through the BH-bulge mass re- 
lation of McLure & Dunlop (2002), which we assume to be 



invariant with redshift (cf. Ho 2007). The completeness of 
the BH mass function depends on both the flux limit of SDSS 
and the widths of the broad emission lines. The adopted BH 
mass function of Vestergaard et al. (2008), complete only for 
M. ;> 1O 8 2 M , yields = -22.4 mag, which corresponds to 
Lq = 10 45 5 erg s" 1 if bright quasars have an average Eddington 
ratio of - 0.2 (Shen et al. 2008). 




z 

FIG. 1 . — Evolution of (a) the luminosity density, (b) quasar BH mass den- 
sity, and (c) duty cycle. In calculating U, the contribution from the infrared 
band, which represents reprocessed emission, has been subtracted. The solid 
lines represent the best linear-squares fits used for calculating r\. In panel b, 
the red square (pf° = 10 1 3 M Mpc- 3 ) is the mass density of local active BHs 
(Greene & Ho 2007). The final results (rf), within the errors, do not depend on 
the particular analytical fits of p% (z) and 8(z). Our calculations assume the 
following Poisson statistical uncertainties: A^q/^q = 0.1, A^q/^g = 0.1, 
and Api so / = 0.1, except for the red point in (b), for which we adopt 

A P r/pr=o.2. 



Table 1 Duty Cycles 



at 


Ml =-22.3 


M* = -22.4 


M* = -22.5 


a 


-4.10 


-4.01 


-3.92 


a\ 


3.34 


3.32 


3.29 


<*l 


-1.60 


-1.59 


-1.57 


^3 


0.26 


0.26 


0.25 



Figure 1 shows the results for U(z), p1 h °(z) and 5(z). Panel 
(b) gives the best fit for /o2 s °(z)- Since the duty cycle S(z) de- 
pends on M^, we plot the results for three choices of (the 
dependence on Lq is absorbed into M R ). We fit the variation of 

S with z using a polynomial log S = Y^=o a i z '> tne coefficients a,- 
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FIG. 2. — Evolution of (a) the luminosity density, (b) quasar BH mass density, and (c) duty cycle. In calculating U, the contribution from the infrared band, 
which represents reprocessed emission, has been subtracted. The solid lines represent the best linear-squares fits used for calculating r\. In panel b, the red square 
(pf = 10' 3 M Q Mpc- 3 ) is the mass density of local active BHs (Greene & Ho 2007). The final results (rf), within the errors, do not depend on the particular 
analytical fits of pt°(z) and <5(z). Our calculations assume the following Poisson statistical uncertainties: A^'q/'I'q =0.1, A^'g/'J'g =0.1, and Api so / =0.1, 
except for the red point in (b), for which we adopt Ap^°/p^° = 0.2. 



are given in Table 1 . We find that 5(z) is very similar for the 
samples of Dahlen et al. (2005) and Ilbert et al. (2005), indicat- 
ing the robustness of the present results. The large uncertainties 
in the photometric redshifts for low-z galaxies preclude us from 
deriving 5 directly for z < 0.3; in this regime, S is extrapolated 
from the analytical fits at higher redshift. We hope to remedy 
this situation with future surveys. 

Our main results are summarized in Figure 2a. We find (1) 
rj ss 0.3 at z ~ 2, and (2) a strong cosmological evolution, with 
r\ decreasing to ^0.03 by z ~ 0. The low value of 77 in the local 
Universe agrees with the results from the continuity equation 
of the BH number density (Merloni & Heinz 2008). The high 
values of 77 at large z are consistent with the estimate from the 
cosmic X-ray background (Elvis et al. 2002), which is mainly 
contributed by AGNs at z « 1 — 2. At this epoch, the spin an- 
gular momentum of BHs may originate mainly from the orbital 
momentum in major BH coalescences (i.e. BH binaries with 
mass ratio ~ 1; those with mass ratio < 1 are referred to as 
minor coalescences) (Hughes & Blandford 2003). 

The strong evolution in rj provides important constraints on 
the accretion history of BHs. There are two ways for BHs to 
grow (Berti & Volonteri 2008): (1) major and minor BH coales- 
cences and (2) accretion of gas. While major coalescences do 
contribute some angular momentum to BHs (Hughes & Bland- 
ford 2003), accretion is the main contributor to their masses 
(King et al. 2008). An initially non-rotating BH rotates very 
fast after accreting about one-third of its mass with the same 
sign of angular momentum (Thorne 1974), and its radiative ef- 
ficiency approaches the maximum value of ?i max = 0.42. How- 
ever, the net angular momentum contributed to accreting BHs 
depends on whether the accretion is prograde or retrograde, as 
described below. 

4. BH GROWTH BY RANDOM ACCRETION 

Supposing that accretion occurs randomly and episodically 
in quasars, then, at episode i with angular momentum A£, and 
mass Am,, we have (as in a random walk), after n episodes, the 
net angular momentum AL. = J2l=i an d = YTi=] ^1 = 



nA£ 2 , whereas the growth in mass AM. = Y^Li = "Am. 
The specific angular momentum of an accreting BH, 



a oc 



y/EA£ 



oc n 



-1/2 



(7) 



nAm 

decreases with the number of accretion events. Consistent with 
the numerical simulations of King et al. (2008), we expect a 
to decrease with decreasing z- The evolution of 77 offers strong 
observational support for this picture. The episode number n 
can be estimated by 

«=-r— = ^)-r—, (8) 
Am Afrj 

where Afo is the episode lifetime and tc(z) is the age of the 
Universe at redshift z- For z ~ 2, we have 5 » 10~ L5 (Fig. lc), 
?c = 2.0 Gyr, and Afo =1.0 Myr as estimated from the prox- 
imity effect of z ~ 2 quasars (Kirkman & Tytler 2008), and 
hence n w 10 2 . The parameter a decreases by a factor ~ 10 
from z ~ 2 to z « 0. On the other hand, random accretion may 
occur through disks whose masses are limited by self-gravity 
(King et al. 2008) to AM sg w (H/R) M„ where H is the height 
of the disk at radius R. We have 



n = f 



R 



10 , 



(9) 



where H/RmlO~ 2 and / = AM. /M. w 1 is the factor by which 
the BH mass has been increased due to accretion. The rough 
consistency between these two estimates of n suggests that the 
random accretion events originate from a self-gravitating disk. 
Additionally, the gas accreted onto BHs is only a tiny fraction 
of the gas content of the entire galaxy (~ 10~ 3 ), so the net an- 
gular momentum from accretion should be independent of the 
total angular momentum. The evolution of the radiative effi- 
ciency found here does not favor accretion with consistently 
prograde angular momentum, as would occur if the typical ac- 
cretion event has much larger aligned angular momentum than 
the BH (Volonteri et al. 2007; Berti & Volonteri 2008). Rather, 
this study strongly favors a picture in which most of the mass 
in BHs was gained through random accretion. 

Accretion rates of ~ 10- 100M© yr" 1 are needed to gen- 
erate the prodigious luminosities seen in quasars. Such large 
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amounts of fuel cannot be sustained by internal sources (e.g., 
stellar mass loss); they almost certainly require external trig- 
gers, in the form of major gas-rich galaxy mergers and inter- 
actions. How do major mergers lead to episodic, random ac- 
cretion? Numerical simulations and direct observations (see 
Schweizer 1998 for review) show that gas-rich mergers drive 
large quantities of gas to the circumnuclear region of the merger 
remnant, where a large portion of the fuel gets consumed in a 
starburst. The mechanisms by which gas gets further funneled 
to smaller scales relevant for feeding AGNs are unclear (see 
Wada 2004 for a review). But two things are certain: the gas in 
the circumnuclear region is clumpy, and in order for the gas to 
accrete, it must get rid of most of its angular momentum. This 
can be accomplished by clump-clump collisions, especially in 
a turbulent environment. Such a process is inherently random, 
and the material that manages to fall in carries little or no mem- 
ory of its original angular momentum. Thus, despite the fact 
that the original source of gas derives from a major merger, the 
final, individual accretion events are essentially random, small 
fragments. 

5. SUMMARY AND CONCLUDING REMARKS 

We derive a simple formalism to link the radiative efficiency 
parameter r\ with the luminosity density, BH mass density, and 
duty cycle of quasars. Applying this formalism to existing 
survey data of quasars and galaxies, we show that rj evolves 
strongly with redshift, with values of <~ 0.3 at z ~ 2 and ~ 0.03 
at z ~ 0. A straightforward interpretation of this trend is that 
the BHs in the quasar population begin with high spins at early 



times, but they gradually spin down to the present day. We sug- 
gest that accretion events onto BHs are inherently random and 
episodic, and that this mode of accretion naturally reduces their 
spin. 

We end with a couple of corollary remarks. 

From the expression for r](z), we can derive p,(z) = 
U(z)/r)(z)c 2 , the cosmic history of BH accretion. Figure 2b 
shows that p.(z) closely traces psfr(z), the cosmic history of 
the star formation rate density. In agreement with earlier stud- 
ies (e.g., Richstone et al. 1998), this suggests that the history of 
BH growth generally follows the history of stellar mass build- 
up in galaxies, at least for redshifts up to ~ 2. 

The physical mechanism responsible for the generation of 
strong radio jets in radio-loud AGNs remains largely a mystery. 
BH spin is often implicated as a critical parameter (e.g., Sikora 
et al. 2007). The principal result of this study is that the average 
spin of BHs in quasars decreases with decreasing redshift. If a 
large spin is required to produce strong jets, we would naively 
expect the fraction of radio-loud quasars to rise toward higher 
redshift. This is the opposite of what is actually observed in the 
SDSS sample of quasars studied by Jiang et al. (2007). 
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AGN group. The research is supported by NSFC via NSFC- 
10325313 and 10733010, and by CAS via 10521001, KJCX2- 
YW-T03, and the 973 project. 



REFERENCES 



Berti, E., & Volonteri, M. 2008, ApJ, 684, 822 
Bongiorno, A., et al. 2007, A&A, 472, 443 
Chokshi, A., & Turner, E. 1992, MNRAS, 259, 421 
Croton, D., et al. 2006, MNRAS, 365, 11 
Dahlen, T., et al. 2005, ApJ, 631, 126 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 430, 604 
Elvis, M., Risaliti, G., & Zamorani, G. 2002, ApJ, 565, L75 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Gebhardt, K., et al. 2000, ApJ, 539, L13 
Greene, J. E., & Ho, L. C. 2007, ApJ, 667, 131 
Ho, L. C. 2007, ApJ, 669, 821 

Hopkins, A. M., & Beacom, J. F. 2007, ApJ, 651, 142 

Hopkins, P. F, Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 731 

Hughes, S. C, & Blandford, R. D. 2003, ApJ, 585, L101 

Ilbert, O., et al. 2005, A&A, 439, 863 

Jiang, L., Fan, X., Ivezic, Z., Richards, G. T., Schneider, D. R, Strauss, M. A., 

& Kelly, B. C. 2007, ApJ, 656, 680 
Kelly, B., Vestergaard, M., & Fan, X. 2009, ApJ, 692, 1388 
King, A. R., & Pringle, J. E. 2006, MNRAS, 373, L90 
King, A. R., Pringle, J. E., & Hofmann, J. A. 2008, MNRAS, 385, 1621 
Kirkman, D., & Tytler, D. 2008, MNRAS, 391, 1457 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Magorrian, J., et al. 1998, AJ, 115, 2285 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K, Maiolino, R., & Salvati, M. 

2004, MNRAS, 351, 169 
McLure, R. J., & Dunlop, J. S. 2002, MNRAS, 331, 795 



Merloni, A., & Heinz, S. 2008, MNRAS, 388, 1011 
Richards, G. T., et al. 2006, AJ, 131, 2766 
Richstone, D., et al. 1998, Nature, 395, 14 

Schweizer, F. 1998, in Saas-Fee Advanced Course 26, Galaxies: Interactions 
and Induced Star Formation, ed. R. C. Kennicutt Jr., et al. (Berlin: Springer- 
Verlag), 105 

Shen, Y, Greene, J. E., Strauss, M. A., Richards, G. T., & Schneider, D. P. 2008, 
ApJ, 689, 160 

Sikora, M., Stawarz, L., & Lasota, J.-P. 2007, ApJ, 658, 815 

Silk, J., & Rees, M. J. 1998, A&A, 331, LI 

Small, T. A., & Blandford, R. D. 1992, MNRAS, 259, 725 

Soltan, A. 1982, MNRAS, 200, 115 

Thome, K. S. 1974, ApJ, 191, 507 

Tremaine, S., et al. 2002, ApJ, 574, 740 

Vestergaard, M., Fan, X., Tremonti, C. A., Osmer, P. S., & Richards, G. T. 2008, 
ApJ, 674, LI 

Volonteri, M., Sikora, M., & Lasota, J.-P. 2007, ApJ, 667, 704 

Wada, K. 2004, in Carnegie Observatories Astrophysics Series, Vol. 1: 

Coevolution of Black Holes and Galaxies, ed. L. C. Ho (Cambridge: 

Cambridge Univ. Press), 186 
Wang, J.-M., Chen, Y.-M., Yan, C. S. & Hu, C. 2008, ApJ, 673, L9 
Wang, J.-M., Chen, Y.-M., Yan, C. S., Hu, C, & Bian, W.-H. 2007, ApJ, 661, 

L43 

Wang, J.-M., Chen, Y.-M., & Zhang, F. 2006, ApJ, 647, L17 

Wolf, C, et al. 2003, A&A, 401, 73 

Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 



